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ABSTRACT 

SMSS J031300.36—670839.3 (hereafter SM0313—6708) is a sub-giant halo 
star, with no detectable Fe lines and large overabundances of C and Mg relative 
to Ca. We obtained VLT-UVES spectra extending to 3060 A showing strong OH 
A-X band lines enabling an oxygen abundance to be derived. The OH A-X band 
lines in SM0313—6708 are much stronger than the CH C-X band lines. Spec¬ 
trum synthesis hts indicate an [0/C] ratio of 0.02 ± 0.175. Our high S/N UVES 
data also enabled us to lower the Fe abundance limit to [Fe/H]( 3 D),NLTE < —7.52 
(3 a). These data support our previous suggestion that the star formed from the 
iron-poor ejecta of a single massive star Population HI supernova. 

Subject headings: stars: Population III. stars; supernovae;general. individual;SMSS 


J031300.36-670839.3 
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1. Introduction 


Keller et al. (2014) described the discovery of SMSS J031300.36—670839.3 (hereafter 


SM0313—6708) and the analysis of the high resolution spectrum taken with the Magellan 
Inamori Kyocera Echelle (MIKE: [Bernstein et ah||2003 ) spectrograph on the 6.5 m Magellan 
Telescope at Las Campanas Observatory in Chile. The spectrum was remarkable for the 
complete absence of lines of iron or any other metal, apart from Mg I and Ca II. The fact 
that no iron line is evident in the signal-to-noise (S/N) 100 (at A = 4000 A) spectrum 

indicated that the iron abundance in SM0313—6708 was at most [Fe/H] = —7. 2j^ making 
this star at least 40 times more iron-poor than HE 1327-2326 ([Fe/H] = —5.6; [Frebel e~ 


al.| |2005D . 


The only elements whose abundances could be measured were lithium, carbon, 
magnesium, and calcium with H(Li) = 12 -|- log(Li/H) = 0.71 ± 0.10, [C/H] = —2.44 ± 0.10, 
[Mg/H] = —4.30 ± 0.10, and [Ca/H] = —7.2 ± 0.10, respectively. The level of lithium is 
presumably related to the decline from the primordial Spite plateau level of H(Li) = 2.2 
during early post main sequence stellar evolution (Lind et al. 2009) and this evidence, 
together with a = 5125 K and \ogg = 2.3 derived from spectrophotometry, is consistent 
with a star ascending the giant branch for the hrst time. The C, Mg, and Ca abundances 
can be attributed to a low energy Pop HI SNe and measurement of the oxygen abundance, 
in particular, offers an important test of this scenario. To measure the oxygen abundance 
in ultra-metal poor K stars, however, it is necessary to work in the ultra-violet and observe 
the OH lines that persist well after the [01] 6300.3 H line has become undetectable. 


^Based on observations obtained with European Southern Observatory (ESO) telescopes 
(proposal 092.D-0742) 

^logarithm of abundance ratio relative to present day solar surface abundances, [A/B] = 


log(nA/nB) - log(nA/nB)0 
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2. Observations 


In order to measure the OH (and CH) lines between 3100 A and 3200 A, and obtain 
better detection limits on N, Na, V, Co, Ni, Fe, and Cu, we used UVES (Ultraviolet and 


Visual Echelle Spectrograph; Dekker et ah 2000) on the VLT (ESO Very Large Telescope) 
at Paranal in Chile. Between September 2013 and September 2014, 31 exposures totalling 
29.26 hours of on-source integration were obtained. The UVES UV spectrum extended 
from 3060 A to 3860 A and the two red spectra from 4786 A to 6803 A. 


The data were reduced using the default settings of the UVES pipeline as implemented 
in the “reflex” reduction environment. There was a high order sinusoidal-like modulation of 
the continuum remaining in much of the UV spectrum. The modulation with an amplitude 
of about 20 % and frequency about 50 A was removed by hand, leaving a residual variation 
of amplitude of 2 % — 4 % and spacing of about 1.5 A that makes the measurement of shallow 
features somewhat uncertain. It was very clear, however, that there were many strong OH 
lines present between the observational cutoff around 3060 A and 3200 A, and that the OH 
lines were much stronger than the CH lines, except at the band head of the CH C-X band 
near 3144 A. In Fig. [^we compare the UV spectra of SM0313—6708 ([Fe/H] < —7.2) and 
HE 0107—5240 ([Fe/H] = —5.4), a C-rich ultra-EMP star with almost identical effective 
temperature and gravity, and possibly at the same evolutionary stage. The total absence 
of metal lines in SM0313—6708 (except for the Mg I lines around 3830 A) is very marked 
as is the different distribution of the forest of molecular lines below 3250 A. The CH lines 
dominate the spectrum of HE 0107—5240 whereas OH lines dominate in SM0313—6708. 
The complete absence of absorption lines in SM0313—6708 between 3250 A and 3800 A 
enables more stringent limits to be placed on those species with strong lines in that region. 






3100 3200 3300 3400 3500 3600 3700 3800 


Wavelength 

Fig. 1.— The normalized UVES UV spectra of SM0313—6708 and HE 0107—5240. Note 
the different distribution of the OH and CH lines in the two stars and the complete absence 
of metal lines in SM0313—6708 except for the similar strength Mg I lines. 
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3. Analysis 


We have carried out a spectrum synthesis of the CH, OH, and NH molecular lines 


in SM0313—6708 using Turbospectrum (Alvarez & Plez 1998 Plez 2012). We adopted 


an interpolated “standard” metallicity MARCS ID LTE model atmosphere structure 


(Gustafsson et al. 2008), for = 5125K, log^f = 2.30, [M/H] = —5.0, but with actual 
abundances from this paper for the atomic and molecular equilibria calculations. We used 


the new and ^'^CH line list from Masseron et al. (2014), while for OH and NH we used 


extensive line lists recently recomputed by Masseron (in prep.). The Masseron OH gf values 


are very similar to those of Gillis et al. (2001) but the NH gf values are signihcantly smaller 


by 0.3 — 0.7dex than the commonly used values of Kurucz (2011). The solar composition 


was adopted from Asplund et al. (2009). 


We rehtted the CH A-X and B-X lines in the blue MIKE spectrum and determined a 
carbon abundance of [C/H] = —2.41 ± 0.05, very similar to the previously determined A-X 
value. The CH C-X band synthesis is in good agreement with the A-X and B-X bands. The 
synthesised spectrum well reproduces the observed CH lines in strength and position for all 
J values and for all branches observed. 

The OH lines extend from 3060 A to redward of 3200 A. Because of the lower S/N 
and the uncertain continuum placement at the shortest wavelengths where the lines are 
strongest, we derived the oxygen abundance in SM0313—6708 by htting the weaker OH 
lines between 3178 A and 3200 A. Fig. I shows clearly that the oxygen abundance is well 
constrained at [0/H] = —2.00 ± 0.10. In HE 0107-5240, the similar temperature and 


gravity star with abundances [M/H] —5.4, [C/H] = —1.58 (Christlieb et al. 2004) and 

[0/H] = —3.0 ( Bessell, Christlieb fc Gustafsson| 2004), these longer wavelength regions 
could not be used for OH htting as they are dominated by GH lines as seen in Fig. We 
note that the treatment of scattering (as done in Turbospectrum), is very important for the 
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analysis of the UV lines in SM0313-6708 and metal-poor giants in general. 

Fig. shows the comparison between the observed and synthetic spectra at the shorter 
wavelengths in SM0313—6708 and HE 0107—5240. The panels on the left show the region 
that was used to measure the oxygen abundance in HE 0107-5240, and the panels on the 
right shows the region around the CH C-X band head. Note the significantly different 
relative strengths of OH and CH lines in the two stars. 

There is no indication of any NH features centered on 3360 A, but the higher S/N of 
the UVES spectrum enables us to reduce the upper limit to [N/H] < —4.2, at least two 
orders of magnitude lower than C and O. 


3.1. 3D—ID corrections to C and O abundances 


The carbon, nitrogen and oxygen abundances reported above were derived from ID 
LTE analysis of the CH and OH lines. One of us (RC) generated a 3D time-dependent 
hydrodynamic model atmosphere of a red giant star with stellar parameters very similar 


to the ones of SM0313—6708 using a custom version of the Stagger code (Nordlund & 


Galsgaard 1995 Nordlund et ah 2009) as well as a ID model atmosphere corresponding to 


the same parameters and relying on the same micro-physics and opacity data. Synthetic 
LTE CH and OH spectra for a range of C and O abundances were produced using the 


SCATE code (Hayek et al. 2011) with 3D and ID models, accounting for the effects of 


radiative scattering; the best fits to the observed spectrum were determined in both cases 
and the 3D—ID abundance corrections were then derived. These 3D—ID corrections 
are: —0.14dex for the carbon abundance as derived from the CH band and —0.53dex 
for oxygen from OH (the negative sign meaning the 3D-derived abundances are lower). 
The 3D-corrected [0/C] ratio in SM0313—6708 is therefore -1-0.02 with an uncertainty of 
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Fig. 2.— Region between 3178 A and 3196 A in SM0313—6708. Observations (thick black 
line); synthetic spectra for [C/H] = —2.41 and oxygen abundance [0/H] = —2.10 (red line) 
and higher abundance [0/H] = —1.90 (blue line). Although the OH lines here are weaker 
than the 3060 A — 3200 A region, the S/N is higher and the continuum level is better defined. 
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Wavelength 



Wavelength 


Fig. 3.— A spectral region showing mainly OH lines (left) and a region aronnd the CH 
C-X band head (right). Top: SM0313—6708. Observations (thick black line); synthetic 
spectra for abundance [0/H] = —2.10 (red line) and higher abundance [0/H] = —1.90 (blue 
line). [C/H] = —2.41 in both syntheses. The 3110 A — 3118 A region has low S/N and 
uncertain continuum placement. Bottom: HE 0107—5240. Observations (thick black line); 
left: synthetic spectra for best htting abundances [C/H] = —1.20, [0/H] = —2.55 (red line) 
and higher abundance [0/H] = —2.40 (blue line). Right: [0/H] = —2.55, [C/H] = —1.20 
(red line) and lower abundance [C/H] = —1.40 (blue line). 
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0.175 dex. A 3D—ID correction of —0.50 dex to the ID nitrogen abnndance was also derived 
with similar 3D and ID LTE calcnlations bnt based pnrely on the comparison of synthetic 
3D and ID spectra. 


The poorer fit of the deep OH and CH lines in Fig. is likely dne to limitations with 
the adopted ID model atmosphere. Onr compntations of CH, NH, and OH lines have been 
carried out in LTE, but in low-density atmospheres it is possible that non-LTE processes 


may be important for cooling the outer stellar layers (Lambert et ah 2013) and for affecting 
the strength of molecular bands. As no study of non-LTE CH, NH, and OH line formation 
is available in the literature, however, we are unable to quantify non-LTE effects on CNO 
abundances derived from these molecular bands. 


3.2. 12^/130 


In order to put limits on the ratio we used two regions of the spectrum where 

the ^^CH and ^^CH lines are well separated and free from blending. There were five strong 
lines in the B-X region and nine in the A-X region. Sections of spectrum centered on each 
of the lines were extracted and summed in both the observed spectrum and the synthetic 
spectra. The summed observed and synthetic ^^CH lines were in good agreement but no 
^^CH feature was evident in the summed spectrum. This allows us to place a three sigma 
lower limit of about 40 for the ratio. This low value of ^^C is consistent with 


massive Population HI star yields with no or little primary nitrogen production (Heger & 


Woosley 2010). 
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3.3. Revised upper limits on the abundances of Fe, Na, Si, Sc, V, Co, Ni, Cu 

Two of the Fel transitions with the largest log 5 / valnes in the observed spectral range 
are located within the highest S/N portion of the UV spectrnm, at 3719.94 A and 3737.13 A. 
Sections aronnd the position of these lines were extracted, snmmed and compared with the 
snmmed synthetic spectra compnted for a range of Fe abnndances. Again, no featnre was 
evident in the observed spectrnm and we estimate a revised abnndance limit for [Fe/H] of 
abont —7.5. 

The non-detection of lines of Sc, Ti, V, Co, Ni, and Cn in the UVES spectrnm in the 
range from 3100 A to 3700 A provided revised lower npper limits for many elements given 
in Table (Frebel et ah, in prep.) The non-detection of the Si I line at 3905.52 A in the 
bine MIKE spectrnm indicates an npper limit of —5.6dex. 

The UVES red spectrnm was in good agreement with the MIKE spectrnm for the 
Mg I lines and the Li I line, bnt the better S/N and weaker tellnric lines in the UVES red 
spectrnm indicates that the Na limit conld be lowered to abont —5.7 dex. Table [T] lists the 
cnrrent best estimates for the chemical abnndances of SM0313—6708. 


Discussion 


Heger & Woosley (2010) have followed the evolution and explosion of metal-free stars 


with masses 10 Mq to 100 Mq and determined their nucleosynthetic yields, light curves, and 
remnant masses. Such stars would have been the first to form after the Big Bang and could 
have left a distinctive imprint on the composition of the early universe. The evolution of 
these stars may be affected in several ways by their zero-metallicity primordial composition. 
Because these stars suffered little mass-loss during their pre-explosive lifetimes, their masses 
at their deaths were larger than more metal-rich stars. The models produces some carbon. 
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but for masses less than about 40 Mq primary production of nitrogen is suppressed. These 
stars ended their lives as compact blue supergiants making their envelope more tightly 
bound, with two interesting consequences: i) The more compact envelope increases the 
amount of mass that falls back after the explosion, meaning that they were more likely 
to collapse to black holes at lower masses compared to modern stars, trapping most of 
the heavy elements; ii) mixing by the Rayleigh-Taylor instabilities is suppressed in these 


compact stars, so it is mainly unmixed material that escapes in the explosion. |Heger 
& Woosley ( 2010| htted their predicted nucleosynthetic yields to several examples of 
metal-poor stars and found that the best fits were obtained from relatively low energy 


explosions of 0.6 — 1.8 B (1 B = 10^^ erg) as compared to hypernova models (Nakamura et ah 


2001). Overall, the large sample of metal-poor stars in Cayrel et ah (2004) was well fit with 


a normal Salpeter IMF including all stars from 10 Mq to 100 Mq (but with a preference for 
stars lighter than 15 Mq), a small amount of fallback and with reduced mixing of 1 % of the 
helium core mass (for reference, the default mixing value adopted to solar metallicity stars 


is 10% of the helium core mass; Heger & Woosley 2010). For the two ultra-iron-poor, but 


high 0/Fe, stars HE 0107—5240 (Christlieb et ah 2004; Bessell, Christlieb & Gustafsson 


2004) and HE 1327—2326 (Aoki et ah 2006 Frebel et ah 2006), the best fits were for similar 


low energy explosions, a lower upper mass range (10 Mq to less than 30 Mq), small fallback 
and an even lower mixing of 0.025. 


The abundances of SM0313—6708 are much more extreme than the Cayrel et ah (2004) 


EMP stars or the two C-rich ultra-EMP stars. The biggest difference are the large carbon, 
oxygen and magnesium overabundances relative to calcium, and the non-detection of all 
other metals. Of particular significance is the oxygen to carbon ratio in SM0313—6708 of 
[O/CjsD = -1-0.02 compared to [O/CjsD ~ —1.0 (Collet et ah 2006) in HE 0107—5240. 


Fitting the current SM0313—6708 abundances (Table 1) against models from Heger 
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Table 1: Chemical abundances of SM0313—6708. 


Element [X/H]id,lte 

[X/H]id,NLTE [X/H](3d),LTE 

[X/H](3d),NLTE [X/H]3d,LTE 

Li I* ** 0.70 ±0.10 

0.73 ±0.10 

0.56 ±0.10 0.71 ±0.10 


C (CH) " 

-2.41 ±0.05 




-2.55 ±0.09 

N (NH) ” 

<-4.2 




<-4.70 

0 (OH) ” 

- 2.00 ± 0.10 




-2.53 ±0.15 

Nal ” 

<-5.7 

<-5.65 

<-5.79 

<-5.66 


Mg I ” 

-4.45 ±0.07 

-4.03 ±0.03 

-4.53 ±0.09 

-4.08 ±0.03 


All ” 

<-6.2 





Si I ” 

<-5.6 





Call 

-7.46 ±0.10 

-7.27 ±0.10 

-7.60 ±0.10 

-7.26 ±0.10 


Sell" 

<-5.5 





Till" 

<-7.1 





VII" 

<-4.3 





CrI 

<-6.3 





Mnl 

<-5.8 





Fel " 

<-7.80 

<-7.28 

<-8.27 

<-7.52 


Col " 

<-5.4 





Nil 

<-6.9 





* 

* 

1 —1 

O 

<-5.5 





Znl 

<-3.4 





Sr II 

<-6.9 





Ball 

<-6.1 





EuII 

<-2.9 






*The abundance of Li is expressed as A(Li) = \og{N(Li)/N{H)) + 12 

** Abundances derived in this paper. 

Estimated (Icr) observational uncertainties in the quoted upper limit abundances are 0.3 dex. For Fe the 
uncertainty is 3cr. Bold print indicates our recommended values used in the fitting of explosive nucleosynthesis 









& Woosley (2010) using their Starfit algorithm]^ finds good fits for progenitor masses 
in the range 40 Mq — 60Mq and energies of 1.5B — 1.8B with low amounts of mixing, 


< 2%) mixing. This low amount of mixing is consistent with Heger & Woosley (2010) 


and hydrodynamic models of Pop III SNe by Joggerst, Woosley & Heger (2009). The best 
matching mass is sensitive to the upper limits of the Na and A1 abundances. Generally, 
lower values on the upper limits to the abundances would push the fit toward higher initial 
masses, stars that produce a stronger odd-even effect. The low limits on N in SM0313—6708 
constrain the mass on the upper end. 


The range of [0/C] abundances produced in a range of core collapse supernovae is 
shown in Fig. [^together with a 40 Mq fit from Starfit. Significantly, the observed [0/C] 


ratio is lower than what Pop III pair-instability SNe of Heger & Woosley (2002) can 
produce or what we can produce for high explosion energies (i7 > 3 B) in more massive stars 
(M > 20 Mq). The structure and pre-SN evolution of stars with mass above 15 Mq is very 


sensitive to initial mass (e.g., Sukhbold & Woosley 2014) and these stars may experience 
chaotic nuclear burning. As a result, the fraction of the C and O layers that are ejected 
varies strongly for a given explosion energy and is rather sensitive to mass and explosion 
energy. Currently there is no good first principle calculations providing reliable predictions 
for explosion energies for a given stellar structure for masses above 15 Mq. 


Our new 3D LTE values provide definite constraints on the [0/C] ratio instead of only 
an upper limit. The new lower value of [Fe/H] derived here strengthens the constraint not 
to eject iron group elements from the inner core of the progenitor star. This limits the 
explosion energy and the magnitude of mixing due to Raleigh-Taylor instabilities during 
the SN explosion. 


^see also http://starfit.orgi (Chan et ah 


2015) 
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Fig. 4.— Top: Range of [0/C] abundances produced in a range of core collapse supernovae 
for three different explosion energies, 1.2 B, 1.8 B, and 3B and mixing of 1.5% (of helium 
core size; see 


Heger & Woosley 2010). The dark gray line and the light gray shading show 


the derived SM0313—6708 and its la error bar of 0.175 dex. Red dots show the [0/C] ratios 


for Pop III pair-instability supernovae (Heger & Woosley 2002). Bottom: The yields from 


stars in the 40 Mq — 60 Mq range provide the best £t to the observed abundance pattern of 
all the elements (see the main text). The 40 Mq fit is shown here. 
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Ishigaki et al. (2014) have proposed an alternative scenario to explain the nnique 


abnndances in SM0313—6708. They investigated snpernova yields of metal-free (Popnlation 
III) stars and find that the high [C/Ca] and [C/Mg] ratios and npper limits of other 
elemental abnndances are well reprodnced with the yields of core-collapse snpernovae (which 
have normal kinetic energies of explosion of IB) and hypernovae (10B) of Popnlation III 
25 Mq or 40 Mq stars. Their best-£t models assnme that the explosions undergo extensive 
matter mixing and fallback, leaving behind a black hole remnant. In these models, Ca 
is produced by static/explosive O burning and incomplete Si burning in the Population 
III supernova/hypernova, in contrast to our suggestion that the Ca originates from the 
hot-CNO cycle during pre-supernova evolution, which is, remarkably, consistent with the 
present data. 


Although Ishigaki et al. (2014) well fits the observed abundances, their modelling 


uses several independent free parameters. For example, their models require signihcantly 
more mixing than is consistent with hydrodynamical simulations of Population III stars 


(Joggerst, Woosley & Heger 2009) and the fallback is adjusted independently of the 


explosion energy. Jetted or asymmetric explosions may be a possibility, however, this is 
another free parameter. In the (Heger fc Woosl^2010) models the mixing is consistent 


with hydrodynamical simulations of Joggerst, Woosley & Heger (2009) and the fallback is 


computed self-consistently with the explosion energy. Fuller discussion of the abundances 
of SM0313-6708 and SN modelling will be presented in a future paper. 


5. Conclusions 

The ultraviolet spectrum of SM0313—6708 is dominated by the lines of the OH A-X 
band which are noticeably stronger than the CH C-X lines. The new line lists of Masseron 
fit the OH and CH lines very well in both wavelength and strength. We could detect 
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no lines of ^^CH, which places a limit of more than 40 on and no lines of NH 

putting nitrogen at least two orders of magnitude less than C and O. This supports there 
having been no mixing from the H-burning shell to the surface within SM0313—6708 and 
there being little nitrogen relative to C and O in the material from which SM0313—6708 
formed. The clear lack of any metal lines apart from Mg I in the UVES ultraviolet spectrum 
also enabled us to lower the abundance limits for those elements whose strong lines are 
normally seen in this region, including Fel, for which we now derive an [Fe/H] upper limit 
of —7.52dex (3 a) from {3D) LTE analysis. We derived an [0/C] ratio of +0.02 with a 
l(j uncertainty of 0.175 dex for SM0313—6708, updated abundances measurements for Mg 
and Ca, and derived new upper limits for nitrogen, the Fe-group elements and some of the 
alpha-elements. These abundances remain consistent with a single massive Population III 
progenitor in the mass range 40 Mq — 60 Mq, low mixing, and modest explosion energies; 
the Ca production remains consistent with an origin from hydrostatic hydrogen burning. 
Improved upper limits on intermediate-mass elements, in particular Na and Al, will further 
help constrain the models in the future. The strong Fe II lines in the vacuum UV between 
2340 A and 2625 A should allow one to constrain - or detect ~ the iron abundance of 
SM0313—6708 down to as low as [Fe/H] = —9. 
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